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Abstract 



We study the one-dimensional Cahn-Hilliard equation with an additional 
driving term representing, say, the effect of gravity. We find that the driving 
field E has an asymmetric effect on the solution for a single stationary domain 
wall (or 'kink'), the direction of the field determining whether the analytic 
solutions found by Leung [J. Stat. Phys. 61, 345 (1990)] are unique. The 
dynamics of a kink-antikink pair ('bubble') is then studied. The behaviour 
of a bubble is dependent on the relative sizes of a characteristic length scale 
E , where E is the driving field, and the separation, L, of the interfaces. For 
EL 3> 1 the velocities of the interfaces are negligible, while in the opposite 
limit a travelling- wave solution is found with a velocity v oc E/L. For this 
latter case (EL <1) a set of reduced equations, describing the evolution of 
the domain lengths, is obtained for a system with a large number of interfaces, 
and implies a characteristic length scale growing as (Et) 1 ^ 2 . Numerical results 
for the domain-size distribution and structure factor confirm this behaviour, 
and show that the system exhibits dynamical scaling from very early times. 
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I. INTRODUCTION 



When a system is quenched from a homogeneous high-temperature phase into a two- 
phase region, domains of the new equilibrium phases form and evolve with time. The 
dynamics of systems described by a conserved order parameter, such as binary alloys or bi- 
nary liquids undergoing phase separation, are conventionally modelled by the Cahn-Hilliard 
equation [I]. At late times after the quench, the domain coarsening is well described by a 
scaling phenomenology, with a single characteristic length scale L(t). 

Recently there has been much interest in the dynamics of phase separation in the pres- 
ence of an external driving field 0-§J, as this has applications to, for example, spinodal 
decomposition in a gravitational field. 

In order to capture the dynamics of a system in the presence of an external driving field, 
an order-parameter-dependent diffusion coefficient (or 'mobility') is required The 
resulting modification of the Cahn-Hilliard equation has been studied by several authors 
both analytically and numerically 0,|5]-|9|]. In all the two-dimensional simulations the two 
phases form structures which align along the direction of the field, with different time- 
dependences for the typical length scales parallel and perpendicular to the field. Analytic 
work on interfacial dynamics has revealed an instability in surfaces perpendicular to the 
field 0,|T31 which explain this configuration. However, given the complexity of the problem, 
no satisfactory scaling results have been obtained. 

In order to obtain some qualitative understanding, we study in this paper the Cahn- 
Hilliard equation in one dimension. This system has been studied in the absence of a field 
by several authors |TT]-|l^] . In the presence of a field, analytic solutions have been obtained 
for a single interface in an infinite system ||. In section |lj we show, however, that these 
solutions are not unique. The direction of the field relative to the kink profile determines 
whether the analytical solution obtained by Leung || is unique, or just one of a family of 
solutions. The work in this section motivates the introduction of an additional characteristic 
length scale inversely proportional to the field strength M. 
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In section [Hj the stationary profile of kink-antikink pairs is studied numerically. The 
behaviour of the of these systems is dependent on the relative values of the characteristic 
length E~ x and L, the distance between the interfaces. These systems are considered in two 
separate limits, EL <C 1 and EL ^> 1. For EL <C 1, the kink-antikink pair form a travelling 
wave, moving at a speed proportional to E/L. For EL ^> 1, the domain walls are essentially 
frozen. 



The limit EL 1 is reconsidered in section [TV], in the context of a system with many 
interfaces. We consider the case in which all the domain lengths, Li, satisfy EL{ <C 1, and 
show that a simple effective dynamics for the interfaces can be constructed, in which the 
velocity of a given interface is proportional to the difference of the inverse lengths of the 
domains on either side of it. This then leads to a simple equation for the time evolution of 
the domain lengths themselves, from which it follows immediately that any characteristic 
length scale must grow as y/t. Numerical simulations of the effective dynamics demonstrate 
scaling, with an average domain size growing as \ft as expected. 

Section |V| concludes with a discussion and summary of the results. 



II. SINGLE-INTERFACE SOLUTIONS OF THE CAHN-HILLIARD EQUATION 

The phase-separation dynamics of a conserved system are modelled by the Cahn-Hilliard 
equation |] 

where E[cp] is the usual 4 Hamiltonian, but with an additional term to include the effects 
of the external field: 

F[<j>] = J d d x Q(V0) 2 + i(l - 2 ) 2 - E.x^j . (2) 

Throughout this paper only the deterministic Cahn-Hilliard equation is considered, i.e., 
the Langevin noise term has been omitted. This restricts the limits of validity of these 
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results to temperature T = or, more generally, to T small enough that the effects of 
thermal fluctuations may be neglected. 

If the mobility A is independent of the order parameter, inserting (Q) into ([!]) shows 
immediately that the term involving E drops out of the dynamics. For non-trivial (and 
physical) results, therefore, the inclusion of an external field requires that the dependence of 
the mobility A on the order-parameter field <fi should be explicitly taken into account P,|5]J§. 
In this paper we take the simplest form for A that maintains the symmetry under — > — 0, 
i.e. A = (1 — a<f) 2 ). Hence rewriting equation ([!]) in one dimension, and defining aE = E, we 
obtain 

d<j> _ | £ d(0 2 ) _ oj 

dt dx 2 dx dx ' 

where 

/! = 0(0 2 -l)-0" 

is the chemical potential and J is the current. Note that we have only kept the leading- 
order term in a, as is conventional. Technically, this is equivalent to taking the limits a — > 0, 
E — > oo, holding E = aE fixed. 

For an infinite system, analytic solutions to equation (H) have been found for a single 
stationary interface ||. The order-parameter profile for a kink solution is given by 

4>{x) = 0oo tanh(0 oo x/v / 2), (4) 

with 

4L = l + V2E, (5) 

valid for both signs of E. The corresponding antikink solution is obtained by inserting an 
overall minus sign in (f|), and replacing E by —E in (||). [This reflects the symmetry of (|3p 
under <p — ► — </>, E — ► —E.] 

In general, solutions to non-linear equations are not unique. The aim of this section is 
to investigate the solutions of equation (|3|) numerically, and to present some simple analysis 
which accounts for the (at first sight surprising) results. 
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The continuum equation (|3|) was discretized on the interval (— L/2, L/2), and antiperiodic 
boundary conditions, <p(—L/2) = —<j)(L/2), were imposed. A kink solution was sought by 
using an antisymmetric tanh profile, with 0(0) = and <p(L/2) > 0, as the initial condition. 
A range of initial values of 0(L/2) was used. This initial state was then evolved under the 
dynamics (|3|) until the system reached a stationary profile. A system size L = 500 was used 
with a mesh size of 0.5. This value of L was found to be large enough that the results showed 
no finite-size effects. The time step used in the iteration was At = 0.005. 

It was found that the number and form of the stationary profiles found was strongly 
dependent on the sign of the field. For each negative field value a unique stationary wall 
profile was found, independent of the initial conditions used, i.e. of the initial value of 0(L/2). 
The final value of 0(L/2) obtained was found to be below the analytically predicted value 
(1 + \/2E)i [equation @]; the discrepancy vanished, however, as the mesh size was decreased 
(see figure 1). 

For a positive field a family of solutions was found, with the final value of 0(L/2) equal 
to the initial value. The interface was found in general to be broader than for E < 0, the 
asymptotic exponential decay to the value 0(L/2) occurring at a rate determined by the 
magnitude of the field, the decay rate decreasing as the field decreases. For larger values of 
the field, this exponential tail became oscillatory (see figure 2). 

Some simple analysis of the differential equation (|3]) can be used to explain some of these 
features. The existence of a solution representing a single interface in an infinite system 
is dependent on the choice of boundary conditions. For a physical solution, two of the 
boundary conditions imposed must be 0(0) = 0, which fixes the position of the interface, 
and 0"(O) = 0, a consequence of the anti-periodicity. The other two boundary conditions, 
0'(O) and 0(oo), must be chosen to ensure the existence of a solution. This places constraints 
on the the possible values of these variables. By the use of some simple stability analysis it 
is possible to deduce the type of constraints this places on 0'(O) and 0(oo) and hence the 
number of values for which a solution exists. 

Equation (H) can be integrated once to give a third order differential equation, 
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Erf? + (30 2 - 1)0' - 0"' = Ecft , (6) 

the integration constant Erf) 2 ^ = Erf) 2 (±00) ensuring that rf>(x) tends to a constant value at 
infinity. Linearizing around this value for large positive x, rf>{x) = 0^ + rf>(x), gives 

0'" + (1-30^)0' -2Erf> oo } = 0. (7) 

It follows that rf> will have solutions of the form exp(Ax), where the possible values of A are 
the roots of the cubic equation 

A 3 + (1 - 30^)A = 2Erf )oo . (8) 

The nature of the roots of this equation is dependent on the values of (poo and E. 

Note that in the present discussion we are taking rf>(oo) > 0, and considering both signs 
of E. This means that we are dealing with a kink, not an antikink. Later, when we consider 
systems with both kinks and antikinks, we will fix £ > 0. Since, however, the equation 
of motion (^) is invariant under E — > —E, rf> — > — 0, it follows that the properties of an 
antikink for E > are identical to those of a kink for E < 0. In the following stability 
analysis, therefore, we restrict the discussion to kinks. 

There are four cases to consider, as illustrated in figure 3. 

a) E > and (Erf)^) 2 > (0^ - |) 3 

Two roots are complex with negative real parts, and the third root is real and positive. 

b) E > and (Erf)^) 2 < (0^ - |) 3 

All roots are real - one positive and two negative. For E — > 0, 0oo — > 1 and the roots 
are A w ±^2, -E. 

c) E < and (Erf)^) 2 < (0^ - i) 3 

All roots are real - one negative and two positive. For E — ► 0, the roots are A ~ ±v2, 
—E. 
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d) E < and {E^f > (0^ - |) 3 

Two roots are complex with positive real parts, and the third root is real and negative. 

For x — > oo, the function <j)(x) has the form 

3 

(f>(x) = / ] dj exp(Ajx). 
i=i 

The coefficients dj are fixed by the initial conditions and the value of 0^. Linear stability 
of the putative solution with <p{x) — *> 0oo as x — > oo requires that vanish at infinity, i.e. 
the coefficients Oj corresponding to eigenvalues \ with positive real part must vanish. 

To see what this implies, consider integrating the third-order equation (^) with the 
initial conditions 0(0) = = 0"(O), which follow from the antisymmetry of the solution, 
and some fixed 0'(O). One also needs to specify the value of the parameter 0^ appearing 
in the equation. The two parameters 0'(O) and 0oo must be varied until the solution 0(x) 
approaches the same value (p^ as x — > oo. If the input values deviate slightly from the 
correct values, the solution will ultimately diverge from 0^ due to the presence of unstable 
modes (with Re A, > 0) in the linearised solution for 0. To ensure that such diverging 
terms are absent, the parameters 0'(O) and 0OO have to be adjusted until the corresponding 
coefficients Oj vanish. This is the familiar 'shooting method' of solving non-linear differential 
equations. 

For E > there is only one unstable mode. Hence for solutions to exist only one 
coefficient must vanish, imposing a constraint of the form ai(0'(O), 0oo) = 0. This defines a 
line in parameter space for which this coefficient vanishes, resulting in a range of possible 
values for 0^. One therefore expects a family of solutions, parametrised by (p^, for E > 0. 

For E < 0, however, there are two roots with positive real parts. Therefore the coefficients 
of both unstable modes are required to vanish. Hence the only choice for the parameters 
0oo and 0'(O) for which a solution exists is at the intersection of the two lines defined by 
setting both coefficients to zero, i.e. ai(0'(O), 0^) = = a 2 (0'(O), 0oo). This results in a 
unique solution for E < 0. 
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The values of E and (poo determine whether the roots of the cubic equation @ are real 
or complex. The existence of complex roots with negative real parts for E sufficiently large 
and positive [case (a) in figure 3] results in oscillatory behaviour. 

For a given stable solution, the eigenvalue A with the least-negative real part determines 
the rate of the asymptotic exponential decay towards 0^. For E small and positive, this 
eigenvalue is —E. Hence for E > 0, the value of E determines a characteristic length in 
the system, E^ 1 , which defines the width of any boundary effect at an interface. This 
characteristic length is important when considering systems with more than one interface 
as in the following sections. For E < 0, by contrast, there is only one negative eigenvalue, 
which approaches a non-zero limiting value for E — > 0: the width of a kink for E < 
therefore remains of order unity for small E. 

How do we interpret the analytic kink solution (£|) in the light of the foregoing discussion? 
For E < 0, the analytic solution is the unique solution we predict. For E > 0, we have argued 
for (and numerically demonstrated) a family of solutions. Which one corresponds to equation 
(f|)? It is easy to show that the exponential decay of (Q) to (p^ as x — > oo is governed by 
the more negative of the two negative eigenvalues arising from the stability analysis. The 
exact solution (f|) therefore corresponds to the case where the amplitude a corresponding to 
the least negative eigenvalue vanishes, i.e. to just one member of the family of solutions. 

Throughout this section we have considered only kink solutions (<p > at positive in- 
finity) and allowed E to take either sign. We conclude by interpreting our results for the 
physically relevant case where the sign of E is fixed (we take E > 0) and both kinks and 
antikinks are present. From the symmetry of the dynamics under — > — 0, E — > —E, it 
follows that for E > there is a family of kink solutions, but a unique antikink solution. In 
a system with well-separated (relative to E~ l ) kinks and antikinks, however, continuity of 
the function will select the kink solution that matches smoothly on to the unique antikink 
solution. When the kink-antikink the separation is small compared to E , more interesting 
behaviour is possible. This is the subject of the next section. 
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III. SINGLE-BUBBLE SOLUTIONS OF THE CAHN-HILLIARD EQUATION 



In this section solutions to the Cahn-Hilliard equation for a kink-antikink pair in a 
periodic system are considered. The kink-antikink pair is a bubble of plus phase moving 
in a sea of minus phase. The velocities of the interfaces were investigated for various field 
strengths and kink-antikink separations. 

The work in the previous section produced a characteristic length E^ 1 , which determines 
(for a kink) the width of any boundary effect present at an interface. In the systems consid- 
ered in this section there is a second important characteristic length, L, the distance between 
the two interfaces. The behaviour of these systems is dependent on the relative values of 
these two characteristic lengths. In this section we will consider two limits, firstly the case 
where EL ^> 1 and then the opposite limit, EL <C 1. The results from the second regime 



motivate a simple model for the many interface problem which is dealt with in section IV. 



A. EL > 1 

The numerical solutions found for a kink-antikink pair in this limit are stationary. The 
value of the order-parameter between interfaces was dependent only on the value of the field 
and was numerically found to be consistent with (f> « ±(1 — \/2E)2, the value given for an 
antikink by the analytical solution (£|). Some examples of the profiles found are shown in 
figure 4. 

These results can be explained using the results of section If we assume that the 
system will reach a steady state and the order-parameter may be written as a travelling 
wave, = 4>(x — vt), then the Cahn-Hilliard equation (^) may be integrated once to obtain 

E<p 2 + (30 2 - 1)4>' - <p"' = J - v<f>, (9) 

where J Q is a constant of integration and v is the wave velocity. In the limit EL ^> 1 
the order-parameter is approximately constant both inside and outside the bubble, ~ 
±(1 — y/2E)% . Inserting this into the above equation in these regions, and demanding that 
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it hold inside domains of both sign, it is clear that one must have Jq = E(l — \p2E) and 
v = 0. 

Therefore the solution is stationary, and since in this limit the separation L of the 
interfaces is large compared to E -1 , the profile is the combination of two single interface 
solutions. 

For one of the interfaces there exists a family of solutions for a fixed value of the field 
whereas the other interface has a unique solution with a plateau height given by ±(1 — \/2E)^. 
This specifies the boundary conditions for the first interface, hence picking out one from the 
family of solutions. 

B. EL -C 1 

In this limit the system was found numerically to evolve to a steady state with a fixed 
profile moving at a finite velocity. This may be shown by a comparison of the order- 
parameter with the current, since for a travelling wave solution equation (|3|) gives —v(jJ = 
—J'. Therefore the order-parameter will be a multiple of the current J, plus an overall 
constant. This expectation is confirmed by the data presented in figure 5. 

As in the previous limit, a long-range boundary effect exists at the kink interface, with 
the order-parameter decaying exponentially to a fixed value outside the bubble. However, 
since EL <C 1 the full exponential tail on the the other side of the interface (i.e. between the 
two interfaces) is missing, and the order-parameter varies roughly linearly between the two 
interfaces, through a height difference proportional to the field. This linear regime between 
the interfaces is evident in figure 5. 

Using the assumption (verified numerically) that the solution may be written as a trav- 
elling wave, equation @ becomes 

E (0 2 - <&) + (30 2 - 1)0' - <T + v (0 + 0oo) = 0. (10) 
where —(poo is the value of the order-parameter in the minus phase far from the interfaces. 
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The integration constant Jo in @ was fixed in (|10D by the requirement 0(±oo) = —(f), 
Integrating this equation over the whole system yields an expression for the velocity: 



where the final result follows from noting that ~ 0^ inside the bubble and — 0oo outside. 

The function (0^ — 2 ) is peaked at the interfaces. This can be used to explain why a 
finite velocity is expected. Consider first the case EL 3> 1, for which we showed above that 
the velocity is zero. For this case, therefore, J dx (0^ — 2 ) must be zero. Since the kink 
and antikink are well separated, and the order parameter saturates at the value 0OO in the 
intermediate plateau region, j dx (0^ — 2 ) is a sum of kink and antikink contributions. 
One of the interfaces, corresponding to the antikink, is known analytically [from equation 
(|)], and the contribution it makes to the integral is 2\/20oo- The contribution from the kink 
interface, therefore, must give an overall negative contribution of — 2y / 20 O o- This is possible 
because, as is clear from figure 4, the kink solution overshoots 0^ and approaches the plateau 
at 0oo from the 'wrong' side. The function 0^ — 2 for the kink therefore contains both 
positive and negative regions, as shown in figure 6. The negative regions evidently contain 
the larger area [note that, although the magnitude of 0oo — 2 in these regions is 0(E), the 
decay constant is also 0(E), leading to an area of 0(1)], in such a way that the integral for 
the kink interface is precisely — 2y / 20 O o- 

Consider now the case EL <C 1. In this case, the interior of the bubble is too small to 
contain the full exponential tail from the kink. The contributions from the kink to j (0^— 2 ) 
will therefore be less negative than when EL 3> 1, resulting in a net positive value for the 
integral and, via (|TTD, a positive bubble velocity, v oc E/L. Therefore in this limit the 
interfaces evolve to a fixed profile which moves with a velocity v oc E/L. The qualitative 
tendency of smaller bubbles to move faster than larger ones was noted by Yeung et al., on 
the basis of numerical studies 0. 

We can make this argument quantitative in the limit of small E, for which 0oo — > 1. 
In this limit, the area under the kink peak in figure 6 approaches the same value, 2y/2, 




(11) 
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as the area under the antikink. The negative region between the two peaks is completely 
suppressed, as the decay length, 1/E is much larger than the separation, L, of the peaks. 
The negative contribution from the semi-infinite region to the left of the kink peak, however, 
retains the value — 2y / 2. The total area under the kink is therefore zero, and the area under 
the kink-antikink pair is 2\/2. Equation ( |TTD then gives a velocity v = \/2E/L for E — > 0. 
Applying this result to the data of figure 5, for which E = 0.01 and L = 20.88 (measured 
between the zeros of 4>) gives v = 6.77 x 10~ 4 , which compares well with the measured 
velocity of 6.82 x 10~ 4 . 

The results of this section will be used in the following section to derive the set of reduced 
equations employed to simulate a multi-domain system. As a first step in generalising this 
approach to many domains, we express the single-bubble calculation in the limit EL <C 1 
in terms of currents. From equation (|]), the current J can be written, up to an overall 
constant, as 

J = —fi' — E<f) 2 = <p'" + (1 - 30V - E<f) 2 . (12) 

The velocity of an interface is related to the discontinuity in the current at the interface 
||, u = AJ/A0, where A J and A0 are the discontinuities in the current and the order 
parameter respectively. For a bubble that moves without change of shape (a travelling 
wave), J must have different values J in and J out inside and outside the bubble. To lowest 
order in E, we can write <ft 2 = (j) 2 ^ = 1 both inside and outside, so that ([12]) becomes 



J = <$'" — 2<p' — Em both regions. Outside the bubble, far from the interfaces, (f) is constant, 
so Jout = —E. Inside, the order parameter changes by an amount of order E over a distance 
L, so we estimate 0' = —E/L (and <fi'" is negligible), giving Jj n ~ — E + 0(E/L) correct 
to 0(E). Finally, both interfaces have the same velocity v ~ E/L, in agreement with our 
earlier result ( [TT]) based on the travelling wave solution. This second approach is, however, 
readily generalisable to the multi-domain situation in the small- E limit. 
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IV. MANY INTERFACE MODEL 



In this section a set of reduced equations is derived using the work of the previous 
section to the describe the dynamics of a large number of interfaces in a periodic system. 
This reduced dynamics is numerically simulated, and the average domain length, domain-size 
distribution function and structure factor are calculated. 

We consider the case in which ELi <C 1 for all domains i, where the Lj are the domain 
lengths. From the arguments of the previous section, we expect that, in this limit, <$' ~ 
—E/Li in domain i, and that the current Jj in domain % is given, up to an additive constant, 
by Ji oc E/Li. The proportionality constant in this relation is the same for all domains (an 
extension of the argument given above for the kink-antikink pair [FJ] gives its value as 4\/2). 

Let us define interface i to have domain length Lj to its right, and Lj_i to its left. If the 
interface is a kink (A0 > 0) its velocity t>j is 

^ocsfl-J-V (13) 



For an antikink (A0 < 0) the velocity vf k is given by the same expression, but with a factor 
(-1). If we arbitrarily assign even i to domains of positive 0, the equation for the time 
evolution of the domain length, dLi/dt = {v i+ i — t>j), becomes 



dt \Li-i L{+i / 

where an overall constant has been absorbed into the timescale. It follows immediately that, 
if L is a typical domain length, dL/dt ~ E/L, i.e. L(t) ~ (Et) 1 ^ 2 . 

The dynamics of the system were studied numerically using this set of reduced equations. 
The system was initially prepared with 100, 000 domains with approximately equal amounts 
of plus and minus phase. The field E was absorbed into the time scale, i.e. we set E = 1 in 
(|TJ). The time step used in the iteration was At oc t. This was chosen because we expect 
£5 (i.e. power-law) growth. Then the At oc t leads to AL oc L, i.e. domains typically change 
length by a fixed (small) fraction of the mean length in each update. Similarly, domains 
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were annihilated if the distance between interfaces became smaller than a specified (small) 
fraction of the average domain size. In the calculation of the average domain size and the 
domain-size distribution function, results were averaged over 10 runs. It was found that the 
average domain size, (L(t)), grew as \/t as expected, (figure 7). The simulation results show 
that the distribution function P(L, t) (the fraction of domains which have size L at time t) 
scales well from quite early times (figure 8). The scaling function was found to be linear at 
the origin, with a gaussian tail. However, we were unable to calculate it analytically. 

Figure 9 shows the structure factor. This was taken as an average of 500 runs. The 
scaling collapse is good except near the peak, where the data are noisy. A ln-ln plot (figure 
10) shows that in the limit k(L) ^> 1 the structure factor possesses the expected Porod tail, 
S(k,t) oc k~ 2 . The constant of proportionality is determined by the average domain-wall 
density 0, S(k,t) = 4pk~ 2 . Hence we would expect a 'Porod plot' of k 2 {L)S(k,t) to tend 
to the value 4 at large k(L). The data, although noisy in this regime, are consistent with 
this expectation. 

The small-fc data in figure 10 suggest the quadratic dependence, S(k,t) oc k 2 . We can 
show that this is the correct small-fc behaviour, following the methods used to establish the 
fc 4 small- k behaviour for the non-driven Cahn-Hilliard equation in dimensions d > 1. 

For k — > 0, equation ([D can be written in Fourier space as 

^ = iEktf) k , (15) 

where the term of order k 2 is negligible and has been omitted. Integrating with respect to 
t, multiplying by and averaging gives the structure factor 

S(k,t) = -S(k,0) + 2(Mt)<P~k(0)) + E 2 k 2 [ dh [ rft 2 ((0 2 [ti]) fc (0 2 [t 2 ])- fc ) . (16) 

Jo Jo 

Since S(k, t) has the scaling form S(k, t) = Lg(kL), where L is a shorthand for (L), it follows 
that S(k, 0) is negligible compared to S(k,t) at late times (large L) and can be neglected. 
Similarly, (0fc(t)0_^(O)) has the scaling form |1J L x h(kL), where A < d/2 = 1/2 follows from 
the Cauchy-Schwarz inequality flBl , so this term can also be dropped for large L. 
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Next we consider the equal-time correlation function D(r,t) = ([0^, — 2 (x)][0^ o — 
4> 2 (x + r)]). It's Fourier transform is the equal-time version of the quantity required in 
(|T6|). The function p(x) = (p 2 ^ — <fi 2 (x) vanishes except near domain walls, so it is essentially 
a domain- wall density function, with delta-function contributions from the walls |15[]. It 
follows that (p) ~ the mean domain- wall density, while D(r,t) has the scaling form 
D(r,t) = L~ 2 f D {r/L). The two-time generalisation is D(r, t 1 ,t 2 ) = L^ 2 f D (r / L x , L 2 /Li), 
where L\ = L{t\). The spatial Fourier transform required in (|16|) has the form T(k, t±, t 2 ) = 
L^ 1 go(kLi, L 2 /Li), which reduces to L\ x h{L 2 jL\) in the limit A; — ^ [we do not expect T 
to vanish in this limit, because p is not a conserved quantity]. 

If L(t) grows as a power of t, the double time integral in ( |l6l) gives, up to constants, 
t 2 /L, so the right-hand side of fllTf ) is of order E 2 k 2 t 2 /L. Requiring that S(k,t) = Lg(kL) 
have the same small-/c form gives S(k,t) ~ k 2 L 3 , and hence L ~ (Et) 1 ^ 2 . It is reassuring 
that this form for L(t) agrees precisely with that obtained from the reduced dynamics (|14|). 

V. CONCLUSIONS 

We have studied, analytically and numerically, the effect of an external driving field on 
the coarsening dynamics of the one-dimensional Cahn-Hilliard equation at T = 0. 

For a single stationary interface, it was shown that the direction of the field for a domain 
wall of a given sign determines whether there is a unique solution or a family of solutions. 
In the latter case, the approach of the interface profile function to its asymptotic value 
is governed by an exponential tail with a decay constant which vanishes linearly with the 
driving field E. There is therefore a new characteristic length scale, E" 1 , in this system. 

The behaviour of a kink-antikink pair ('bubble') falls into two classes, characterised by 
the relative values of the interface separation, L, and the new characteristic length, E~ l . In 
the limit EL 3> 1, the bubble profile is stationary. In the opposite limit, EL «l,a bubble 
of plus phase moves through the minus phase with a velocity v oc E/L. 

For the many-domain coarsening dynamics an equation of motion for the domain lengths 
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was derived [equation ([14])], in the which the length of a given domain changes at a rate 
determined by the lengths of the domains on either side. The mean domain size grows as 
\f~E~i. Despite the apparent simplicity of this model, we have so far been unable to make 
further analytical progress. Numerical simulations, however, demonstrate dynamic scaling 
(for the domain-size distribution and the structure factor), and confirm the predicted \ft 
growth law. An argument for the observed k 2 behaviour of the structure factor at small k 
has been presented. 

In a related work, the 1-D Ising model with Kawasaki spin-flip dynamics, biased in 
one direction, has been studied []17| |. Numerical results were obtained for a range of volume 
fractions, including the case of equal volume fractions simulated here, and y/t growth demon- 
strated. There, as here, rather general arguments lead to \fi growth for all volume fractions 



- for the present model from dimensional analysis of equation ([14]) . For the Ising model, 
however, exact results for the domain-size distribution (and for the 'persistence exponent' 
fT7|j ) were obtained in the limit where one phase occupies a negligible volume fraction, so it 
is interesting to consider the present model in this same limit. 

When one phase (the plus phase, say) occupies a very small volume fraction, the system 
consists of domains of plus phase separated by (typically) much larger domains of minus 
phase. To a first approximation, therefore, each plus domain can be treated as an isolated 



bubble in the sense of section |TJ. Each such bubble then moves at a speed E/L, where L 
is its length (we are taking EL <C 1 here), all bubbles moving in the same direction. As 
a result, small bubbles catch up with larger bubbles, with which they then merge, and the 
combined domain slows down. Just before the bubbles merge, the approximation of treating 
them as independent breaks down, when the size of the intervening minus domain becomes 
comparable with the size of the bubbles. This presumably leads to negligible corrections 
to the scaling functions, however, in the small volume- fraction limit. This new model is 
so simple that analytic progress may be possible. However, we will leave this question for 
future work. 
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FIGURES 



FIG1. Stationary single-interface solutions of equation (|3|). There is a unique solution 
for each E < 0. (a) E = -0.1, (b) E = -0.3, (c) E = -0.5. 

FIG2. Stationary single-interface solutions of equation ([3]). (a) E = 0.2, (b) and (c) 
E = 0.3: two examples from a family of solutions. For each value of the field there is a 
different decay length. 

FIG3. Graphical solution (schematic) of equation @ for the four cases discussed in the 
text. 

FIG4. Kink-antikink solutions with L « 107 and L w 161, for E = 0.1, corresponding 
to the limit EL > 1. 

FIG5. Comparison between the order-parameter profile and the current for a kink- 
antikink pair of size L fa 21, for E = 0.01, corresponding to the limit EL <C 1. 

FIG6. Plot (schematic) of 0^ — (j) 2 (x) for a kink-antikink pair in the limit EL 1. 
The total area under the curve is zero. 

FIG7. Time dependence of the square of the average domain length for the dynamics 
described by equation flli]) . The mean domain size grows as t 1 ^ 2 . 

FIG8. Scaled domain-size distribution function for the dynamics described by equation 
(0). The data correspond to times 10 (□), 20 (+), 40 (*), 80 (O), 160 (x), 320 (A), 640 
(0). 

FIG9. Scaled structure factor for the dynamics described by equation (|14|). The data 
correspond to times between 10 and 640, as in figure 8. 

FIG10. Double-logarithmic plot of the scaled structure factor. Times between 10 and 
640 are shown, as in figure 8. The straight lines in the small- Ac (L) and large- A; (L) regimes 
have gradients 2 and —2 respectively. The latter is the expected Porod regime; an argument 
for the former is given in the text. 
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